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We present an entirely perturbative QCD determination of the crit- 
ical phase transition flavor number N^ r of QCD. The results obtained 
are compared with various determinations of Ny by non-pertrubative 
methods, including lattice QCD. The wider physics implication of the 
existence of the Banks-Zaks regime of QCD with only weakly inter- 
acting quarks, is discussed briefly. 

Keywords: Banks-Zaks infrared fixed point in QCD 
PACS: 12.38.-t, 
E-mail: ndili@uh.edu 




F. N. Ndili 



Physics Department, University of Houston, 
Houston, TX.77204-5506, USA. 



August, 2005 



Abstract 



1 



1 INTRODUCTION 



Gies and Jaeckel in a recent paper [1] raised the question of reliable 
determination of the critical flavor number Nj r separating perturbative and 
non-perturbative regimes of QCD from the high flavor end. They analyzed 
the chiral phase structure of QCD using a functional renormalization group 
method, and obtained a value of NJ = 10.0 ±0.4. A number of earlier deter- 
minations of the same quantity by other authors exist [2-6], based on similar 
considerations of the chiral phase transition of QCD. Non-perturbative lat- 
tice QCD methods have also been used [7] to find Nf . The importance of 
NJ as a critical phase transition parameter of QCD has generally been em- 
phasized in these determinations. By noting that the take off point common 
to all the above formulations and determinations of Nj r is the Banks-Zaks [8] 
observation of a weakly coupled two loop infrared fixed point at high flavor 
number in QCD, and that this feature of QCD can be exploited not only 
non-perturbatively but also perturbatively, we arrive at the possibility that 
NJ can also be determined by entirely perturbative methods. Such purely 
perturbative analysis of the Banks-Zaks (BZ) infrared phase of QCD were al- 
ready considered in another context of searching for the existence of infrared 
fixed points in higher order QCD beta function [9-11]. Our focus in this 
paper is on the original two- loop Banks-Zaks infrared fixed point, presented 
as a tool for perturbative determination of the phase transition critical flavor 
number NJ of QCD. The results so obtained can be compared with var- 
ious determinations of NJ by non-perturbative phase transition methods, 
including lattice QCD. With largely consistent firm values emerging from 
these different sources, we discuss the wider physics implication of the exis- 
tence of the BZ regime of QCD and its only weakly interacting quarks. 



2 The Banks-Zaks Domain of QCD 

A basic observation made by Banks and Zaks in their paper [8] is that along 
side the ultraviolet asymptotically free fixed point of QCD that holds for 
Nf < 33/2, there exists close to this limit point in flavor space, an infrared 
fixed point that is induced on QCD by its two loop perturbative beta function. 
The QCD beta function is given by: 

(3(a) = —ba 2 {l + ca + c 2 a 2 + c 3 a 3 + c 4 a 4 + ■•■—> oo) (1) 

where a is the QCD coupling constant in the form a = a s /7t. The first 
few expansion coefficients b, c, C2,cs are known, and have the specific values 
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[12-15]: 



b = 



33 - 2N f 
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It is seen that the two loop beta function (3(a) = —ba 2 (l + ca) = has an 
infrared fixed point solution at a = a* = — 1/c. The fixed point value depends 
on flavor, varying from a* = 0.012 at N f = 16 to a* = 0.21 at N f = 8.05. 
For Nf < 8, the infrared fixed point vanishes. Because of the smallness of 
the QCD coupling constant in this fixed point domain, Banks and Zaks sug- 
gested that the domain can be exploited by making perturbative expansions 
in the flavor parameter (33/2 — Nf). 

The BZ domain has however other features which others have noted. It 
is chirally symmetric, and contains massless fermions and gluons [1-4]. Be- 
cause of the weak coupling state of the domain, quarks are not confined 
there. The domain has no mass gap or condensate such as would arise from 
chiral symmetry breaking, and it has no physical particle spectrum such as 
form from qq bound states, or Nambu-Goldstone bosons. Finally the domain 
being massless is scale invariant and conformally symmetric [6]. All these 
properties provide a variety of avenues by which one can explore the nature 
and physics of the BZ domain of QCD, and how this domain relates to the 
rest QCD. 

A central point is that if QCD has the BZ domain controlled as it is 
by an infrared fixed point that exists only for certain high values of flavor 
number, there must exist in QCD not only a phase structure, but a critical 
flavor number Nf at which transitions occur between the different phases of 
QCD. Banks and Zaks [8] sketched a phase structure for QCD, but a more 
authentic phase structure was later constructed by Appelquist et. al. [2,3]. 
Based on this QCD phase structure and various non-perturbative methods, 
determinations have been made of Nj r as follows: 

The value of critical flavor obtained by Appelquist et. al. [2,3] from 
their chiral phase structure and non-perturbative gap equation for QCD, is 



Nf = 11.9. 
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Gies and Jaeckel [1] combining a model of chiral phase structure of QCD 
and Nambu-Jona-Lasinio (NJL) type four fermion interactions, find Nf = 
10.0 ±0.4. 

Sannino and Schechter [4] exploring the same chiral phase structure of 
QCD, use a non-perturbative effective potential and perturbative anomalous 
dimensions, to obtain a value of Nf = 11.7. 

Harada, et. al. [5] from analysis of meson masses in large flavor Bethe- 
Salpeter equation obtain a value for Nf = 12. 

Miransky and Yamawaki [6] from a general perspective of conformal phase 
transition for non-perturbative dynamics of gauge fields, including QCD and 
its BZ domain, obtain a value for the QCD critical flavor number Nf = 11.9. 

Iwasaki et. al. [7] in a lattice modelling of QCD chiral phase structure 
and confinement, obtain a value Nf > 7. 

We argue that comparable values for Nf can be obtained from an en- 
tirely perturbative analysis, based on the original perturbative expansions 
suggested by Banks and Zaks in powers of (33/2 — Nf). First we note that 
the methods listed above rely on a suitable modelling or choice of QCD 
phase structure or phase diagram, before one can set up a corresponding 
non-perturbative QCD calculation or simulation for Nf. In the perturba- 
tive approach discussed below, we use only the firmly established fact that a 
perturbative infrared fixed point exists in QCD at two-loop order, and for a 
range of high flavors. 



3 Purely Perturbative Searches for NJ 

The BZ idea is to exploit the two-loop infrared fixed point near Nf < 33/2 
by making perturbative expansions of physical quantities, around the point 
Nf = 33/2 in flavor space. For this, we rewrite the two-loop infrared fixed 
point coupling: 

°* = 4=153fW (33/2 " jV,) <5> 
in the form a* = e(33/2 - Nf) where e = -4/(153 - 197V/. We next eval- 
uate e at Nf — 33/2, obtaining e = 8/321. Finally we define a suitable BZ 
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perturbative expansion parameter as the quantity: 



(8) 



«o = a* mm = — (33/2 -N f ) (6) 

This quantity a Q can have a range of values from a = at Nf = 33/2 to 
a Q = 0.21 at Nf = 8.05 where the two-loop infrared fixed point vanishes. In 
terms of a , the two-loop infrared fixed point eqn(6) becomes reparametrized 
at any one value of a Q , by the power series: 

* 1 a o i 2 i 2 3 i 3 4 i /i-r\ 

a =-c = l-(19/4)a =a0 + Uao + Uao + Ua ° + (7) 
where u = 19/4. Correspondingly : 

(a*) 2 = al+2ual + 3u 2 a 4 o +4u 3 al + 5u 4 al+6u 5 a 7 o + 7u e a 8 o + 8u 7 a 9 o + 9u 8 al + . 
and 

(a*) 3 = a3 + 3Ma 4 + 6M 2 a 5 + 10M 3 a 6 + 15M 4 a 7 + 2lM 5 a 8 + 28M 6 a 9 + 36M 7 a 10 + 

(9) 

Any physical quantity expanded perturbatively in powers of a Q for any 
a D value within the above infrared fixed point range in the BZ flavor space, 
should be a convergent series. While the point Nf = 8.05 provides a lower 
limit below which no such convergence of a perturbative series is expected, 
it can happen that the perturbative series in a breaks down even before the 
point Nf = 8.05 is reached, due to onset of QCD phase transition bound- 
aries, or to specific non-perturbative QCD dynamics. Each point of break 
down of a physical perturbation series in a Q defines a critical flavor number 
Nf r > 8.05 in this NLO case. Thus by studying various perturbative QCD 
series reformulated as expansions in a , and noting the points of break down 
of each series, one can deduce corresponding values of Nf r from equation (6). 
We now consider some standard QCD perturbation series, and find their ra- 
dius of convergence in the BZ flavor space. 

(i) As our first perturbation series we take the ratio quantity: 

Jwwi- ^r^^fT* < io > 

Its QCD perturbation expansion is written: 

1=1 
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where 



R = a + r 2 a 2 + r 3 a 3 + r 4 a 4 + .... (12) 

with [16-20] 

r 2 = 1.986 - 0.115^/ 

r 3 = -6.637 - 1.200JV/ - 0.005iV 2 - 1.2395(EQ / ) 2 /(3EQ 2 ,) 

We evaluate this QCD perturbative series R e + e - in the specific BZ region 
of high flavors and NLO infrared fixed point. For this we set a = a* in 
equation (12); also we rewrite each of the coefficients r\ in terms of a Q using 
equation (6). We obtain 

n = 0.0885 + 4.6144a (13) 



r 2 = -27.798 + 54.77a - 8.05a 2 (14) 

The BZ perturbation series for R e + e - is now obtained by substituting equa- 
tions (7-9) and equations (13-14) into equation (12). We obtain: 

R e+e _ = R( a *) = a Q + 4.84a 2 + 0.2172a 3 - 184.35a 4 - 2131.445a° 

- 17870. 74a^ - 13.01468 x 10 4 a^ - 87.33858 x 10 4 a^ 

- 55.5123 x 10 5 a^ - 33.9406 x 10 6 a£° • • • (15) 

One can next use our requirement that the perturbation series equation 
(15) defined as it is in a fixed point domain, must be a convergent series, to 
place a bound on a Q and correspondingly on Nf. Before discussing equation 
(15) further, we write down the analogous BZ expansions for other QCD 
physical perturbation series. 

(ii) Our second example is the Bjorken Sum rule for polarized electron 
nucleon deep inelastic scattering, whose QCD perturbation series is [21] 

R e Bj = a + k ± a 2 + k 2 a 3 + (16) 

where 

k 1 = — --N f = -0.917 + 13.375a D (17) 

and 

k 2 = 41.4399 - 7.6073^/ + 0.17747iVj = -35.7641 + 70.2495a o + 285.729a 2 . 

(18) 
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By plugging equations (7-9) and (17-18) into equation (16), we obtain 
the BZ expansion for the Bjorken sum rule : 

R e B . = R e ( a *) = a + 3.833a^ - 8.5405^ - 267.2234a 4 - 2533.49^ 

- 18, 929. 7a 6 Q - 12.6894 x 10 4 a^ - 79.8579 x 10 4 a^ 

- 48.1868 x 10 5 a^ - 28.2121 x 10 6 a^° • • • (19) 

(iii) As our third QCD perturbation series, we take the Bjorken Sum rule 
for neutrino scattering on nucleons. Its QCD perturbation series is [22] 

R [ b) =a + ha 2 + f 2 a 3 + (20) 

where 

f x = 5.75 - O.UN f = 1.58 + 17.833a G (21) 

and 

f 2 = 54.2298 - 9.4968^/ + 0.2392iVj = -37.34588 + 64.332a + 385.1117a*. 

(22) 

By plugging equations (7-9) and (21-22) into equation (20), we obtain 
the BZ expansion for the Bjorken neutrino sum rule : 

flM = flM (a*) = a Q + 6.33a* + 18.057a^ - 84.315a 4 , - 1360.3677a^ 

- 11742.8a° - 8.429 x 10 4 a^ - 55.2112 x 10 4 a^ 

- 34.2104 x 10 5 a^ - 20.4055 x 10 6 a^° • • • (23) 

(iv) Our fourth case is the QCD perturbation correction to r hadronic 
decay: 

r(r~ — > v T + hadrons) 
= T(r- -> u T e~u e ) ^ 
Its QCD perturbation series is [20,23]: 

R T = a + T X a 2 + r 2 a 3 + (25) 

where 

19 

n = F 3 - —B x = 0.08325 + 15.21494a G (26) 
with F 3 = 1.9852 - 0.1153iV> ; and fa = (2N f - 33)/6. Also we have 

r 2 = -38.44269 + 107.0978a o + 254.2058a 2 , (27) 
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By plugging equations (7-9) and (26-27) into equation (25), we obtain the 
BZ expansion for the r hadronic decay width : 

Rl h) = R T (a*) = a Q + 4.83325a^ + 0.1223a 3 - 183.36a 4 - 1849.236^ 

- 13, 927.0a° - 9.29142 x 10 A a 7 o - 5.7986 x !0 5 a s o 

- 34.647 x 10 5 a^ - 20.087 x 10 6 a* • • • (28) 

(v) The last QCD perturbation series we consider is QCD correction to 
scalar Higgs boson hadronic decay width described by the QCD perturbation 
series [24]: 

R H (a) = 3[1 + ha + h 2 a 2 + h 3 a 3 + ] (29) 

where hi = 5.66667 ; 

h 2 = 35.93996 - 1.358657V/ = 13.52235 + 54.5158a (30) 

and 

h 3 = 164.139-25.771^+0.25897^) = -190.57975+691. 1551a +416.952a^. 

(31) 

By plugging equations (7-9) and (30-31) into equation (29), we obtain 
the BZ expansion for the hadronic Higgs decay width : 

B$ = R H (a*) = 3 + 17.0a o + 121.37a* + 360.754a 3 + 47.7a 4 - 9.4869 x 10 3 a^ 

- 9.972 x 10 4 a^ - 7.742 x 10 5 a^ - 52.975 x !0 5 a 8 o 

- 33.7675 x 10 6 a^ - 20.449 x 10 7 a^° • • • (32) 

4 Analysis of BZ Perturbation Series of 
Physical Quantities 

We now look at equations (15), (19), (23), (28) and (32) and make our 
deductions. First we comment on the structure of these BZ series. We notice 
that the first few terms of each series have signs opposite to the rest of 
the infinite series. Furthermore, simple ratio test of convergence with these 
early terms does not yield a consistent result. Since we are interested in the 
overall convergence of each infinite series, we can drop any finite number of 
terms without affecting the convergence property of each series. Accordingly 
we shall drop the first four terms of each series and impose our convergence 
requirement from the fifth term upwards. As for ratio test of convergence, we 
use the more sensitive Raabe ratio test [25] . This states that for a convergent 
infinite series of terms: u n ,u n+ i , the quantity : 
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n - 1 > P > 1 (33) 

\U n +l J 

where n is chosen to be n > N and N is a fixed positive integer, chosen 
here to be the fifth term of each series. We can then choose in succession n = 
5, 6, ...10 .... in equation (33), obtaining values of P, and deducing the con- 
straint values on a Q . The corresponding constraint values on Nf yielding Nf T 
are obtained from equation (6). These values of NJ obtained for the various 
QCD perturbation series we considered above, are summarized in Table 1. 
We find a fairly consistent set of values for the critical flavor number Ny . 
These values compare very well with the results quoted earlier from different 
non-perturbative QCD methods. One feature of our result is that the value 
of NJ obtained for each physical quantity studied, decreases gradually with 
increasing n along each series. The trend is that for n — > oo, the value of NJ 
converges to 8.05 like a* itself in eqn.(7). In the process we find we reproduce 
all the results of various non-perturbative calculations quoted earlier. Our 
results do not yield the lower lattice QCD limit of Nj r > 7, reported in [7], 
nor any values Ny < 6 suggested in [9,10]. Other independent investiga- 
tions using Pade approximants [26] place the QCD infrared phase transition 
boundary well above Nf > 6 in agreement with our result. 



5 Summary and Conclusions 

It is certain from several independent calculations and sources that the criti- 
cal phase transition flavor number of QCD lies well above Nf > 8, associated 
with the BZ two-loop infrared fixed point in QCD. 

A remarkable fact coming out of the existence of NLO BZ infrared fixed 
point region in QCD, is that quarks of flavor Nf > 8 if they exist at all, 
can only be weakly interacting, unable to form bound states of mesons or 
hadrons. Therefore they must remain unconfined and observable in their 
naked color and fractional charge. While various searches for conventional 
quarks of fractional charge and naked color appear so far to yield no positive 
results [27], the existence of the BZ region of QCD is sufficiently indicative 
that these quarks may still be found in the form of very high flavor quarks, 
probably in cosmic rays. 

Additionally, if we remain in the infrared low energy region of high flavor 
QCD, the BZ analysis shows that the inter-quark coupling qq weakens with 
increasing flavor variety. Thus quark pairs at Nf are more strongly paired 
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than quark pairs of Nf > NJ meaning : 





ON 



f 



< qq >=< (34) 



If this trend is assumed to hold continuously across the phase transition 
boundary at Nf < Nj r , one may infer a regime of ordering of the inter- 
quark couplings: tt < bb < cc < ss < dd ~ uu. In turn this will reflect 
in observable flavor dependence of quark condensates, and on the ratios of 
these condensates. Dosch [28] for example obtained that 

= 0.65 < 1. (35) 

< uu > 

However not enough data presently exist on the values of < cc >, < bb >, 
< tt > condensates, to check the above BZ extrapolation. Data on these 
higher flavor quarkonium states, qq potentials and decay widths will also be 
helpful. 



BZ series 


n = 5 


n = 6 


n = 7 


n = 8 


n= 9 


Re+e- 


N f > 12.51 


N f > 11.78 


Nf > 11.27 


N f > 10.89 


N f > 10.59 


R Bj 


N f > 12.02 


N f > 11.37 


N f > 10.92 


N f > 10.59 


N f > 10.33 


R Bj 


N f > 12.63 


N f > 11.71 


N f > 11.14 


N f > 10.74 


N f > 10.45 




Nf > 12.06 


N f > 11.34 


Nf > 10.87 


Nf > 10.53 


Nf > 10.27 


K H 


Nf > 13.3 


N f > 12.07 


Nf > 11.37 


Nf > 10.90 


Nf > 10.54 



Table caption: Raabe [25] ratios of convergence equation (33), calculated 
for various values of u n /u n+ i and for different QCD physical perturbation 
series. 
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